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We investigate the thermodynamical aspects of the Casimir effect in the case of plane parallel plates made of real metals. 
The thermal corrections to the Casimir force between real metals were recently computed by several authors using different 
approaches based on the Lifshitz formula with diverse results. Both the Drude and plasma models were used to describe a 
real metal. We calculate the entropy density of photons between metallic plates as a function of the surface separation and 
temperature. Some of these approaches are demonstrated to lead to negative values of entropy and to nonzero entropy at 
zero temperature depending on the parameters of the system. The conclusion is that these approaches are in contradiction 
with the third law of thermodynamics and must be rejected. It is shown that the plasma dielectric function in combination 
with the unmodified Lifshitz formula is in perfect agreement with the general principles of thermodynamics. As to the Drude 
dielectric function, the modification of the zero-frequency term of the Lifshitz formula is outlined that not to violate the laws 
of thermodynamics. 
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I. INTRODUCTION 



Recent advances in experimental investigation of the Casimir effect (see Refs. pL0|| and the review |llj| ) have given 
^ ' impetus to extensive theoretical studies of different corrections to the Casimir force. Casimir force arises between 
two closely spaced neutral bodies due to the existence of zero-point electromagnetic fluctuations. It is one of the rare 
macroscopic manifestations of quantum phenomena. For this reason, it received widespread attention. Moreover, 
currently the Casimir effect finds applications in fundamental physics for constraining hypothetical forces predicted 
by different extensions to the standard model |12 lJ] and also in nanotechnology [@|§. 



Originally the Casimir force was computed between two infinitely large plane parallel plates made of ideal metal 
fL5fl . Corrections to this ideal result are caused by the geometrical factors (restricted area of the plates and surface 
(-h i roughness), finite conductivity of the boundary metal and nonzero temperature. Geometrical factors are detailly 
examined in the literature (see, e.g., their study in Refs. Jlf^-p0| and also in [jll]]). Corrections caused by the finite 



conductivity of a metal j2l|-[26[ and by nonzero temperature [ 23, 271, 28 1 , when considered separately, also received much 



attention and wholly satisfactory results were obtained (see also [[L l[| ) . 



For experimental purposes the combined effect of different corrections to the Casimir force was found to be of large 
importance. The effect of surface roughness, combined with any other corrections, can be effectively computed by 
the method of geometrical averaging j pT|,p!9[ . So one comes face to face with the problem how to find the combined 
action of finite conductivity and nonzero temperature onto the Casimir force. At first glance it would seem that there 
is an easy way to solve this problem. Use could be made of the famous Lifshitz formula |pl| , p9[ for the Casimir force 
at nonzero temperature acting between two dielectric semispaces by the substitution of the dielectric permittivity 
function describing real metals (on the basis of the plasma model, Drude model or optical tabulated data for the 
complex refractive index). This was done recently by different authors po|-^8| and unexpectedly led to conflicting 
results. 

In Refs. IH,^,^] the corrections to the Casimir force due to the combined effect of the finite conductivity and 
nonzero temperature were calculated in the framework of the Lifshitz formula and of the free electron plasma model. 
The obtained results are in agreement. Temperature corrections are positive and smoothly transform to those for an 
ideal metal in the limit of infinite plasma frequency. 

In Refs. |30],[u) the Drude dielectric function was substituted into the Lifshitz formula. The temperature corrections 
to the Casimir force were found to be negative at small space separations between the plates. The asymptotic values 
of the Casimir force at high temperature (large separations) were found to be two times smaller than for ideal metal 
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irrespective of how high the conductivity of real metal is. Thus, the results of |MJ,[3l| do not convert smoothly to the 
results, given by the plasma model, when the relaxation goes to zero and also to the results, found for an ideal metal, 
when the plasma frequency goes to infinity. So extraordinary properties of the obtained results were attempted to 
explain in fl30| , [3l| by the principal role of nonzero relaxation. Mathematically these properties are caused by the zero 
value of the reflection coefficient at zero frequency as given by the Lifshitz formula for photons with perpendicular 
polarization in the framework of the Drude model. 

In Refs. |37],|38| both the plasma and the Drude models were used supplemented by the special prescription modifying 
the zero-frequency term of the Lifshitz formula in the same way as was done in Ref. [^2) for the case of an ideal metal. 
As a result, large temperature corrections arise to the Casimir force at small separations that are linear in temperature. 
At large separations (high temperatures), asymptotic of the Casimir force in fl37| , |38[ do not demonstrate any finite 
conductivity correction starting from a separation of several micrometers. 

As explained in Refs. p4p5| ], actually the Drude model is outside of the application range of the Lifshitz formula and, 
specifically, the zero- frequency term of this formula must be modified in an appropriate way in order to incorporate the 
dissipative media. This was demonstrated on the basis of a new derivation of the Lifshitz formula |pd|] in the framework 
of quantum field theory in Matsubara formulation. In Ref. the new prescription for the zero-frequency term of the 
Lifshitz formula was proposed which is not subject to the above-mentioned disadvantages (see also |36|). Discussions 
regarding the correct description of the thermal Casimir force between real metals are, however, being continued (see 
recent Comments f3j],[4(J and Replies J4l],|^j supporting just the opposite points of view). The necessity to conclusively 
resolve a problem is apparent when it is considered that the experiment is already nearing the registration of the 
thermal corrections to the Casimir force. 

In the present paper we analyse the thermal Casimir force between real metals on the basis of the fundamental 
principles of thermodynamics. Entropy for a system of photons between the realistic metallic plates is calculated. It 
is shown that in the approach of Refs. [3C|j3l| ], entropy is negative within the separation range where the temperature 
correction to the Casimir energy density computed in |3(],|3l| is negative. It is proved also that both in Refs. |30 31 
on the one hand, and in ]3^|58| on the other hand the Ncrnst heat theorem, or the third law of thermodynamics, is 
violated. Thus, the approaches of Refs. ||(],|3l],|37],[38| are unacceptable from the thermodynamical point of view. As 
for the results of Refs. J32] , p3| , p5 36| (for plasma model) and Refs. (for Drude model), they are shown to be in 

agreement with the general principles of thermodynamics. 

This paper is organized as follows. In Sec. II the general expression for the entropy of a system of photons between 
metallic plates is presented. Sec. Ill contains the computational results for the entropy in the framework of the Drude 
model using different approaches. In Sec. IV the analogical results obtained in the framework of the plasma model 
are given. In Sec. V the reader finds conclusions and discussion. 



II. ENTROPY FOR PHOTONS BETWEEN PLATES MADE OF REAL METAL 

We consider two plane parallel plates made of real metal which are in thermal equilibrium with a heat reservoir at 
some nonzero temperature T. Let a be the space separation between plates. The modern derivation of the free energy 
per unit area for the system under consideration is based on quantum field theory in the Matsubara formulation and 
^-regularization method. The result is |jLL|| 
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with g ; = (£f /c 2 + fcj^ 2 ) 1 / 2 , fc ; = [e(i(;i){;f / c 2 + k^ 2 ] 1 / 2 are the reflection coefficients for the modes corresponding to 
two different polarizations. Here e is the frequency dependent dielectric permittivity of a plate material computed 
along imaginary frequency axis at discrete Matsubara frequencies £; = iirlksT /h with I = . . . — 2, — 1, 0, 1, 2 . . . , ks 
is the Boltzmann constant, and k± is the modulus of the wave vector component in the plane of the plates. 
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The result ([!]) is obtained by the solution of a one-dimensional scattering problem on the axis perpendicular to the 
plates. In fact, an electromagnetic wave which is coming from the left in one semispace is scattered on the vacuum 
gap between semispaces and there are reflected and transmitted waves (see Refs. for details). 

It is important to keep in mind that the scattering problem leading to Eqs. (|l|), (2|) has the definite solution only 
under the requirement that 

lim CW) = C ^ 0. (3) 

If this is not the case (like for metals described by the Drudc model or for dielectrics, see the next section) some 
additional conditions must be used to fix a solution of the scattering problem at zero frequency. As an example, for 
dielectrics the results ([j]), (||) are restated including the zero- frequency contribution by using the unitarity condition 
and dispersion relation. However, as to the case of the Drude model, describing a medium with dissipation, the 
unitarity condition is not applicable, and, therefore, a solution of the scattering problem at zero frequency remains 
indefinite. Due to this fact, metals described by the Drude model are outside of the application range of formulas (jl]), 
(Q) for the Casimir free energy density at nonzero temperature. The special prescription concerning the zero-frequency 
term of Eq. (|l|) must be introduced in order that the dissipative media could be described on the basis of this equation. 
This prescription must be in accordance with the laws of thermodynamics and other general physical requirements. 

It is easily seen that Eqs. (|l|), (Q) lead to the famous Lifshitz formula |2l],^] for thermal Casimir force between 
two semispaces which is obtained as F(a,T) = ~dFE(a,T)/da. Thus we arrived at the conclusion that the Drudc 
metals are outside of the application range of the Lifshitz formula at nonzero temperature and a special prescription 
is needed in order to describe them in a consistent way. 

In the limit of zero temperature, the Casimir energy density of the zero-point electromagnetic oscillations is reob- 
tained from Eq. (|l|) as 

f- />oc />oo 

E{a) = 4^ J d ^j k±dk± l ln I 1 ~ r " (C ' k±) e ~ 2aq> 

+ ln[l-ri(£,fc ± ) e - 2 ^]}. (4) 

From (Q) the Lifshitz formula for the Casimir force at zero temperature is obtained F(a) — —dE(a)/da. It is a 
nontrivial result that the Lifshitz formula at zero temperature is applicable for nondissipative as well as for dissipative 
media. This was demonstrated in Ref. @| (see also f23f| ) through the consideration of a supplementary electrody- 
namical problem and is explained by the fact that the only point £ = which gives an important contribution to the 
discrete sum (0) does not contribute to the integral (Q). 

Note that the generally accepted terminology "the Casimir energy density and force at zero temperature" is of some 
ambiguity. It is really correct in the sense that the sum of the zero-point energies (not the free energies) is computed. 
This terminology, however, disregards the fact that there is a significant thermal dependence of the energy density 
and force through the thermal dependence of the dielectric permittivity (see the next section). Because of this, when 
one calculates the Casimir energy density at, say, T = 300 K, the value of e(i£, T) is substituted into Eq. (0), not 



e(i^,T = 0) (see, e.g., Refs. 25 2qj3l| , |32| ). Thus, instead of E(a), a more exact notation for the quantity ([4|) would 
be E T (a). 

If we multiply (|l]) or (ji|) by 2irR, where R 3> a is a radius of a sphere at a separation a from the semispace, one 
obtains the Casimir force in the configuration of a sphere near a plate at a temperature T, or at zero temperature, 
respectively Jwl]. 

According to thermodynamics, the entropy per unit area of the system under consideration is 

S(a,T) = ±[E T (a)-F E (a,T)}, (5) 

where Fe(cl,T) is given by Eq. ([I]) and Et(cl) = E(a) from Eq. (|J). In the next sections entropy is calculated for 
the plates made of real metal as described by the Drude or plasma model dielectric functions. In doing so, special 
attention is paid to the zero- frequency term (I = 0) in Eq. (0) and to the fulfilment of Eq. (^|) . 

III. ENTROPY IN THE CASE OF METALLIC PLATES DESCRIBED BY THE DRUDE MODEL 

It is common knowledge that the Drude dielectric function 
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e D (w) = l- * , (6) 

UJ[UJ + tj) 

where oj p is the plasma frequency and 7 <C ui p is the relaxation frequency, gives a good approximation of the dielectric 
properties for some metals, e.g., for aluminum. This approximation was widely used in combination with the Lifshitz 
formula (||) to calculate the finite conductivity corrections to the Casimir force at zero temperature fi|J^ fu]]25|fj6f| . To 
do so the Drude dielectric function along the imaginary frequency axis was considered 

Before we proceed further, we note that the dielectric function (f?J) violates the requirement (||) for any nonzero 
value of the relaxation frequency 7. Due to this, as discussed in the previous section, the Drude metals are outside 
of the application range of the Lifshitz formula (|l ) at nonzero temperature with an unmodified zero- frequency term. 
If, nevertheless, one substitutes Eq. ([?]) into Eq. (lj), as was done, e.g., in Refs. ||6|, |3(i| , |3l| , several questionable results 
follow which are in contradiction with the limiting cases of metal described by the plasma model (see the next section) 
and of an ideal metal (see Introduction and a detailed discussion in Refs. |34|-|3(|). Using the dielectric function of 
Eq. (^), the values of reflection coefficients (||) at zero frequency are the following 

rfj(0,fc_L) = 1, ri(0,fcj_) = 0. (8) 

From the formal point of view, some troubles are connected with the second of Eqs. (||) because for an ideal metal 
ri(0,fcj_) = l. 

In order to present the crucial argument against the substitution of the Drude dielectric function (Q) into the 
unmodified Lifshitz formula (0) at nonzero temperature, we calculate the entropy per unit area given by Eq. (|^) in 
the framework of the Drude model. Thus, as an example, consider Al plates with the parameters |4g] 

uj p » 12.5 eV w 1.9 X 10 16 rad/s, (9) 
7 w 0.063 eV « 9.6 x 10 13 rad/s. 

In Fig. 1 the computational results for the entropy are presented at T = 300 K as a function of the separation between 
the plates. Note that the plasma frequency practically does not depend on temperature. As to the value of the 
relaxation frequency from Eq. (||), it is given for the temperature under consideration. It is seen from the figure 
that the entropy is negative within a wide separation range that it not acceptable from a thermodynamical point of 
view. The separation interval < a < 4.1 /im, where the entropy is negative, coincides with the interval where the 
negative temperature corrections arise in the approach used in f|(],|3l|] (see the detailed discussion in J34|). Negative 
temperature corrections are in conflict with the evident physical arguments (with increase of temperature the number 
of photons in the modes and thereby force modulus should increase) . Here we show that they are also in contradiction 
with the general physical principles. 

As noted in the Introduction, except of the immediate application of the Lifshitz formula in combination with the 
Drude model [ |30"||3l| ] , different prescriptions were proposed in literature modifying the zero- frequency term of this 
formula in the case of real metals. In |5?],[38| it was postulated that 

rjj(0 J fcj.) = ri(0 J fcx) = l, (10) 

as in the case of an ideal metal. This prescription was criticized on physical grounds in |34f| . 

In Ref . the other prescription was proposed which is the generalization of the prescription of Ref. |^2| , formulated 
for an ideal metal. According to [p4| in the framework of the Drude model the reflection coefficients at zero frequency 
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As shown in p3|, prescription ( |TT| ) leads to wholly satisfactory results. 

In order to test all the above approaches for conformity to the general principles of thermodynamics, we find the 
dependence of the entropy (j^) on temperature at some fixed plate separation, say, a = 2 /im. To accomplish this, one 
should take into consideration that except of the evident dependence of Eqs. (pi) and (0) on temperature there is the 
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aforementioned significant thermal dependence of Ed given by Eq. (f7j) through the relaxation parameter 7 = j(T) 
(coinciding with the thermal dependence of resistance). The dependence "f(T) is linear at temperatures higher than 
0.25T*, where T* is the Debye temperature (T* = 428 K for Al [|6|). At lowest temperatures j(T) follows the power 
law T n (n = 2 — 5 depending on the metal). The complete dependence of the nondimensional normalized quantity 
j(T) = 2a"f(T)/c on the temperature is plotted in Fig. 2 by the use of tabulated data for Al (here we neglect the 
small residual resistivity caused by the scattering of electron waves by static defects that is outside of the frameworks 
of the Drude model J47|). 

Now we are in a position to calculate the dependence of entropy on temperature for all the above approaches. 
Calculation was performed using Eq. (||) and also Eqs. (Q), ([|), (0) and Fig. 2. The results are presented in Fig. 3. 
The long-dashed curve is obtained on the basis of the Lifshitz formula with an unmodified zero-frequency term 
(i.e., Eq. (H) was used for the reflection coefficients at zero frequency). Remind that this approach was exploited in 
Refs. pCj , pi| . The short-dashed curve is calculated with the modification of the zero-frequency term of the Lifshitz 
formula in accordance with Eq. ( |io|) (approach of Refs. (3^^]). The solid curve is calculated with the modification 
of the zero-frequency term according to Eq. ( |TT1) suggested in Ref. |34|] . Note that all the above approaches differ by 
the value of the zero-frequency term only for perpendicular polarization. 

As is quite clear from Fig. 3, for both long-dashed and short-dashed curves the values of entropy at zero temperature 



are not equal to zero. In the approach used in Si(0) = —0.5 MeVm 2 K 1 , and in the approach used in |B7 38 

5 2 (0) = 0.016 MeVm^KT 1 . In both cases the value of S(0) depends on the parameters of the system under 
consideration (like the separation between the plates and the plasma frequency) that is in manifest contradiction 



with the third law of thermodynamics (the Nernst heat theorem |t8 49|). It is notable that the entropy density given 
by the long-dashed curve is negative in a wide temperature range (compare with Fig. 1 where the result at a fixed 
temperature is presented). It is easily shown that the values of entropy density at zero temperature, given by the 
dashed curves, are related by 

5 2 ( )-5 1 (0) = ^|, (12) 
Ibna 



where £(3) is the Riemann zeta function. The right-hand side of Eq. (|12| ) is equal to one-half of the coefficient near 
temperature in the zero-frequency term of the Lifshitz formula for an ideal metal. This is because in |5?],[38| the same 
values (|l^) for the reflection coefficients at zero frequency were postulated as for an ideal metal, whereas in p0| , ^l| , 
according to (||), one- half of the result for an ideal metal was used. In fact the correct result for a real metal lies in 
between of these two possibilities. 

Contrary to the dashed curves, for the solid curve (approach of Ref. |54|), 5(0) = and therefore, in accordance 
with the laws of thermodynamics. 

IV. ENTROPY IN THE CASE OF METALLIC PLATES DESCRIBED BY THE PLASMA MODEL 

For the free-electron plasma model the dependence of dielectric function on the frequency is given by 

2 2 
e p (w) = l-^, ep (<0 = l + ^. (13) 

This dependence was widely used to calculate the finite conductivity corrections to the Casimir force at the separations 



of order 1 fim 21-24}]. In Refs. []32 33 it was applied to compute the effect of nonzero temperature and finite 
conductivity in the framework of the Lifshitz formula ([!]). Note that the plasma dielectric function practically does 
not depend on temperature. 

The preference for the plasma model as compared with the Drude one is the fulfilment of condition (||) with the 
dielectric function (|l3|). As a consequence, the scattering problem underlying the Lifshitz theory has a definite solution 
leading to Eqs. (|l|), (^|), (^) including the zero-frequency contribution. The reflection coefficients (|^) at zero frequency 
in the framework of the plasma model take the form 



(ck± - Ju>% + c 2 k\ 
, I • (14) 
ck x + Ju$ + c 2 k 2 X/ 

Note that Eq. (14) can be obtained from Eq. (|ll|) when the relaxation frequency goes to zero. Because of this, in the 



approach of (34 ] the results for the plasma model are obtainable from the results for the Drude model in a limitimg 
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case 7 — » (the results obtained in |3(|[3l| by the use of the Drude model on the basis of Eq. (g) have no smooth 
connection with a plasma model approach; the plasma model by itself is not considered in ]3^,|3l]]). In the alternative 
approach |3^,Q, as distinct from both ]3^-|35[| and [ [30pl| ], the conditions (|l0|) are postulated in the framework of 
both Drude and plasma models. 

Now let us find the dependence of the entropy density on temperature in the framework of the plasma model on the 
basis of different approaches. As an example, Al plates are used once more at a separation a = 2 /im. The value of the 
plasma frequency for Al is given by Eq. (|9|) [ p| . Calculations were performed using Eqs. ( |l|), (p|) , (||), (0) and (|l3|). 
The results are presented in Fig. 4. The dashed curve is obtained in the framework of Refs. ]37| , |3q | using prescription 
@. The solid curve is calculated on the basis of Refs. [|32|-|35| with no modification of the Lifshitz formula, i.e. with 
the zero- frequency reflection coefficients (|l4|). 

As is obvious from Fig. 4, for the dashed curve the value of entropy at zero temperature S(0) = ^(O) = 
0.016 MeVm~ 2 K _1 and is not equal to zero. It depends on the parameters of the system that is in contradic- 
tion with the third law of thermodynamics. In contrast to this, for the solid curve obtained on the basis of the 
fundamental Lifshitz formula, S(0) — as it must be from the third law of thermodynamics. 

In the framework of the plasma model it is not difficult to obtain the analytical expression for the entropy at low 
temperatures fcgT <C ksT e ff = hc/(2a). For this purpose the perturbation expansion of the Casimir energy and free 
energy in powes of two small parameters So /a = c/(aui p ) and T/T e // can be used (see Refs. [|i 33 3j| where these 
expansions arc presented in details; they are applicable for separations a > X p — 2nc/ujp where X p is the effective 
plasma wavelength). Substituting the mentioned perturbation expansions into Eq. one obtains 



S(a,T) 



fcgCg ) ( T 
8na 2 



Jo 



Jeff 

2tt 3 T 
45C(3) T eff 



T 



45C(3) T eff 



(15) 



Here the terms up to (T/T e ff) 3 were included. The powers of So/a higher than one are contained only with powers 
(T/T e ff) n , n > 3. This analytical expression corresponds to the solid curve in Fig. 4. It is seen from Eq. (|l5| ) that the 
entropy approaches zero as the second power of temperature in the framework of the unmodified Lifshitz formula. 
On the basis of the approach proposed in Refs. P4pq] the perturbation expansion of entropy is given by 



g f rji\ fc aC(3) S ( S \ 



(16) 



where S(a, T) is expressed by Eq. (|l5|). The first contribution in the right-hand side of Eq. ( |lrj| ) is the value of entropy 
at zero temperature 



5f fc aC(3) / S \ 



(17) 



At a — 2 /im one obtains from Eq. ( fl7| ) the above value of ^(O). S(a, 0) depends on both the separation between the 
plates a and the penetration depth of the electromagnetic oscillations into the plate material So in contradiction with 
the third law of thermodynamics [p8p9[| . 



V. CONCLUSION AND DISCUSSION 



In the foregoing we have considered the thermodynamical aspects of the Casimir force acting between real metals 
at nonzero temperature. The necessity of considering these aspects stems from the controversial results obtained by 
different authors (see Refs. p0|-p8[) and continuing polemic p9|-p2[. The further importance to this problem is added 
by the rapid progress in experiment on measuring the Casimir force. At the moment, there is a contradiction between 
the experimental results of [Q and the theoretical approach of [3(||3lJ which leads to large negative temperature 
corrections at separation of about 1 ^m (see the discussion in |3^,4^]). On the other hand, the experimental results 
of [§~H do not agree with the computations of |5?],[38| that lead to large (although positive) linear in temperature 
corrections to the Casimir force at separations of about 100 nm (see pijl). 

Bearing in mind that there are many influential factors in so precise experiments, it is highly desirable to offer 
some decisive theoretical arguments providing a way to give preference to one of the theoretical approaches. As 
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shown above, thermodynamics gives the possibility to make a selection and to reject the approaches that are not in 
accordance with the most fundamental physical principles. 

In the present paper we calculated the entropy density for photons between two parallel plates made of real metals 
described by the Drude or plasma models. It is shown that in the approach of Refs. ||(],|3l|, based on a direct 
application of the unmodified Lifshitz formula in the case of Drude metals, entropy takes negative values in a wide 
range of related parameters. The value of entropy density at zero temperature, as given by the approach of p0| , |3l"| | , 
is shown to be nonzero and dependent on the parameters of the system under consideration in contradiction with 
the Nernst heat theorem. These unacceptable properties of entropy confirm the conclusion of Refs. |34|,|35| that the 
dissipative metals described by the complex dielectric permittivity of real frequency are outside of the application 
range of the Lifshitz formula at nonzero temperature. To describe the thermal Casimir force for such metals a special 
prescription should be adopted modifying the zero-frequency term of the Lifshitz formula. 

One prescription of this kind (the same as was proposed in p2| for the case of ideal metal) was suggested in 
Refs. [^,^8|. We show that although the entropy density in the approach of |3^j3^] is positive, the value of entropy 
at zero temperature is not equal to zero and depends on the parameters of the system. Thus, this approach is also in 
contradiction with the third law of thermodynamics. 

One more prescription for the zero- frequency term of the Lifshitz formula was proposed in |34| . It is the general- 



ization of the receipt of 22 for the case of real metals. We show that for the prescription suggested in |34j| , entropy 
is non-negative at all temperatures and takes zero value at zero temperature. Hence the prescription of Ref. |34| is in 
agreement with the general principles of thermodynamics. 

In this paper we report also the results of the computation of the entropy in the framework of the plasma model. 
The plasma model does not take dissipation into account. It belongs to the application range of the Lifshitz formula. 
It is shown that the application of the unmodified Lifshitz formula in combination with the plasma model leads to 
wholly satisfactory results: the entropy density is positive and takes zero value at zero temperature. The application 
of the modified Lifshitz formula, as in Refs. |5?],[58]], leads to the violation of the third law of thermodynamics. 

To conclude, the approach of Refs. |3(],|3l| on the one hand and of Refs. (3^^^] on the other hand must be rejected 
as they are in contradiction with the general principles of thermodynamics. In the case of the Drude metals only the 
approach of Ref. Q| fits thermodynamical requirements. It is also in accordance with the present experimental results. 
As to the case of the plasma metals with no account of dissipation, the unmodified Lifshitz formula is applicable and 
the results of Refs. p2|-p6[ are in agreement between themselves and with the general principles of thermodynamics. 
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FIG. 1. Entropy of photons between Al plates described by the Drude model at T = 300 K as a function of space separation 
computed using the approach of Refs. [30,31]. 
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FIG. 2. Dimensionless relaxation frequency of Al as a function of temperature. 
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FIG. 3. Entropy of photons between Al plates described by the Drude model at a separation a — 2 /im as a function of the 
temperature. Long-dashed curve was computed using the approach of Refs. [30,31], the short-dashed curve was obtained with 
the approach of [37,38], and the solid curve using the approach of [34]. 
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FIG. 4. Entropy of photons between Al plates described by the plasma model at a separation a — 2 /im as a function of the 
temperature. Dashed curve was computed using the approach of [37,38] and the solid curve was obtained with the approach of 
[32-36]. 
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